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5 The enigmatic radio source Sagittarius A* (Sgr A* ) at the centre of our Galaxy appears 
tQ be a low-luminosity version of active galactic nuclei (AGN) in other galaxies^. By analogy 
with AGN models^, it has been proposed that Sgr A* may be a massive accreting black 
hole^'^. This is apparently confirmed by dynamical evidence that the center of the Galaxy 
has a dark object with a mass < IO^^Mq which seems to have a gravitational influence on gas 
and stellar motions out to ~ 1 pc^'^'®. The black hole hypothesis is, however, problematical 
because no model of Sgr A* has been able to explain the observed spectrum in any self- 
consistent way, and there is no consensus on either the mass of the black hole or the mass 
accretion rate®"®. Sgr A* has been observed in the radio^"^^, sub-nim^^'^^, infrared^^"^®, 
and X-ray^®~^® bands, and the various detections and flux upper limits cover more than ten 
decades of frequency, from v < 10® Hz (meter wavelength radio band) up to i^ > 10^® Hz (150 
keV X-rays). We present a robust model of Sgr A* in which a IO^Mq black hole accretes at 
the rate of a few xlO~®M0yr~^. This model fits the entire spectrum self-consistently. The 
unique feature of the model is that the flow is advection-dominated, i.e. most of the energy 
which is viscously dissipated in the differentially rotating flow is carried along with the gas 
and lost through the horizon. The apparent success of this model in explaining the data 
may be considered "proof" that horizons are real and that a massive black hole does exist 
at the Galactic center. 

Figure 1 shows the spectral data on Sgr A* 9~18 -pj^g bulk of the emission occurs in sub-mm at 
v ~ 10^^ Hz and in near infrared at i^ ~ lO^'* — lO^'*^ Hz, but between these two bands there is a 
pronounced and very significant dip in the luminosity in the far IR band^^'^^. (Note that the lO^'*^ Hz 
fiux determination^® is from a crowded field and is better interpreted as an upper limit.) Sgr A* has 
also been detected as a variable soft X-ray source (3-30 keV, variable by a factor of 6 ^^•^'^') and there 
are recent upper limits^® on its luminosity in two hard X-ray bands (35-75 and 75-150 keV). There is 
no direct information in optical and UV because of obscuration^; however, the maximum luminosity in 
these bands is limited to < 3 x lO^^ergs"-'^ .-'^'^'-'^-'^ 

If Sgr A* is to be understood as an accretion powered system, it is necessary to model the viscous 
fiow of the accreting gas as well as the radiation mechanisms whereby the dissipated energy is converted 
into the radiation we observe. The most commonly adopted paradigm in this field is the so-called thin 
accretion disk modeP'® which assumes that the cooling is efficient and that the energy released through 
viscosity is radiated immediately. This leads to a number of simplifications in the dynamics of the fiow 
and in the calculation of the emission. However, the model does not apply to Sgr A* since it predicts 
a spectrum^ which is very different from the data shown in Fig. 1. 

Since the thin disk model makes approximations which may not always be valid, some authors have 
expanded the model to include the effects of radial pressure gradients and radial energy transport. The 
extended theory is referred to as the slim disk model^® and one of its features is that it allows energy 
advection, i.e. the gas is free to transport some of the viscously dissipated energy as stored entropy. 
The fraction / of the energy which is advected is solved for self-consistently as a function of radius 
R. The slim disk equations have been successful in the study of disk instabilities^®'^'^ and in modeling 
boundary layers of accretion disks around white dwarfs and pre-main sequence stars^^'^^. 

The present work grew out of investigations of advectzon-dommated iiows^^~^'^ , where energy advec- 
tion is not a perturbation, or included merely for self-consistency, but actually dominates the physics of 
the fiow. In terms of the parameter /, we consider accretion fiows where / ^ 1, and the local radiative 
efficiency 1-/^1. In the case of accretion onto a black hole, advection-dominated fiows are possible^® 
whenever the mass accretion rate satisfies M < 10~^ — lO~^MEdd, where MEdd = 2.2 x 1{]~^ M MqYI'^ 
is the Eddington accretion rate and M is the mass of the accreting object in solar mass units. The 
model of Sgr A* presented here corresponds to M ~ l{]~'^MEdd- 
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The basic framework of our low-M advection-dominated models is described in ref. 25. Because 
the cooling is inefficient (1 — / ^ 1), these models have a very high ion temperature. The ion pressure 
therefore becomes dynamically important, causing the gas to be almost spherical in morphology and to 
rotate slowly. The model thus represents a hybrid between a so-called "ion torus" ^^ and a non-rotating 
spherical accretion flow^®'®'^, except that we solve self-consistently for the advection of energy and the 
transport of angular momentum. We also compute the thermal structure of the accreting gas more 
completely and obtain realistic spectra. 

We model viscosity by the usual a prescription^^, where the kinematic viscosity coefficient is taken 
to be ac'l/QK, Cs being the local sound speed and Qk the Keplerian angular velocity. This simple 
prescription can plausibly account for the effects of hydrodynamic turbulence, convection^^'^'*, and 
magnetic stresses. We assume that the accreting plasma is composed of gas and magnetic fields, such 
that a fraction (3 of the pressure is supplied by the gas and a fraction 1 — /3 by the fields. We take 
a, j3 to be independent of R. We set a = 0.3 (as in our previous work^®) and discuss two values of 
(3, viz. 0.5 (equipartition between gas and magnetic pressure) and 0.95 (gas pressure dominant). The 
parameters a and (3, along with the mass M and accretion rate M , complete the specification of a 
model. 

Given these parameters, and using a self-similar solution^^ as the local representation of the fiow, 
we solve^® for all the properties of the accreting gas, such as the density p, angular velocity il, radial 
velocity vr, sound speed Cg, etc., at each radius R. We determine the ion and electron temperatures, 
Ti and Tg, self-consistently by solving for energy balance among various processes like viscous heating, 
advection, energy transfer from ions to electrons (via Coulomb collisions), and radiative cooling of 
the electrons. For the cooling we include synchrotron emission, bremsstrahlung, and Comptonization. 
One simplification is that we use a Newtonian model of gravity down to the Schwarzschild radius, 
Rs = ^GM/c^ , and assume that the gas disappears through the black hole horizon at i? = Rs- Also, 
we do not model the details of the sonic transition^® near the horizon or the gravitational redshift of 
the escaping radiation. These approximations will doubtlessly lead to quantitative errors in the results 
but we believe that the model captures most of the esssential physics. 

The solid lines in Fig. 2 show the variations of Ti, Tg and 1 — / as functions of R for one of our 
models, with a = 0.3, /3 = 0.5, M = 10^ Mq and M = 1.3 x W'^MEdd = 2.9 x lO-^Mgyr-i. At large 
radii, Ti and Tg are both very nearly equal to the local virial temperature, Tyir ~ 2 x 10^^ l3(Rs / R) K. 
Energy transfer from ions to electrons is very efficient here and so the two species come into thermal 
equilibrium. However, the cooling of the electrons is inefficient, and the fiow is therefore advection- 
dominated. For R < 10^ Rs, the electron temperature saturates at Tg ~ 10^'^K but the ion temperature 
continues to track Tyir- Although electron cooling is efficient at these temperatures, the ion-electron 
coupling becomes weak, and once again the fiow is advection-dominated. Thus, at all radii, the radiative 
efficiency 1-/^1. 

The solid line in Fig. 1 shows the spectrum predicted by this model. Basically, there are four peaks 
in i>L(i>), where L(i>)di> is the power radiated by the source between frequencies i> and i> + dv. The 
peak labeled S is due to synchrotron emission by the thermal electrons in the magnetic field of the 
plasma, the peaks CI and C2 are due to Comptonization^*^ (single scattering and double scattering) 
of the synchrotron radiation by the hot electrons, and the peak B is due to bremsstrahlung radiation. 
The presence of these four peaks is a robust feature of all low-M advection-dominated models, but the 
frequencies at which the peaks appear and their relative heights depend on the parameters. 

The position of the synchrotron peak S depends primarily on M , while its height varies roughly 
as M . Thus, merely by fitting the sub-mm data and requiring that the model satisfy the far-IR limits 
at i^ ~ 10^^ — lO^'* Hz, it is possible to estimate M and M fairly well. The sharp drop at the high-i^ 
end of the S peak arises because most of the luminosity is emitted at the innermost radius where the 
synchrotron emission is highly self-absorbed. Allowing for fiuctuations in the magnetic field strength 
and gravitational redshift effects, and modeling the spectrum more accurately above the self-absorption 
limit, will round the peak and produce a more gradual fall-off, but a significant dip between S and CI 
is always expected. 

At radio frequencies, v < 10^*^ Hz, the model spectrum lies below the observations. This emission 
comes from large radii where the gas temperature may conceivably be modified by radiative transfer or 
other transport effects which are not included adequately in the model. Another possibility is that some 



of the accreting gas may undergo shocks and participate in an outflow, as advection-dominated flows 
seem prone to do^^'^'*, and this gas may radiate in radio waves through shock-accelerated non-thermal 
electrons. In any case, the luminosity at these frequencies is very low and the discrepancy is probably 
not serious for the overall viability of the model. 

Because the synchrotron emission is highly self-absorbed, the emission at each v in the S peak 
originates essentially at a single radius. The radiation at 86 GHz, for instance, comes from R = blRs = 
1.1 AU. This is in excellent agreement with a recent measurement^"'^ '^^ which gives a size ~ 1.1 AU. 

The CI and C2 peaks are obtained by single and double Compton scattering of the synchrotron 
photons. The frequency shift between the S and CI peaks and that between CI and C2 are determined 
primarily by Tg since the mean energy boost {A) of a photon in each Compton scattering depends on the 
distribution of Lorentz j of the scattering electrons. In the models shown here, the maximum electron 
temperature is Tg^max = 8 — 14 x lO^K which gives (A) ~ 30 — 100. Note that the electron scattering 
optical depth is very small, Tgg < 10~^. Therefore, most of the synchrotron photons escape and only a 
small fraction ~ Tgg is Compton-scattered; the fraction that is scattered twice is ~ t^^. The heights of 
the CI and C2 peaks relative to S are therefore proportional to {A)Tgs and {A)^Tgg respectively. The 
flux in the S peak is ex M, as is the optical depth Tgg. Therefore, the flux in the CI peak varies as 
~ M^ and that in C2 as ~ M^. If the source has a variable M for any reason, these peaks will vary in 
amplitude, with the largest fluctuations expected in C2. Indeed, the luminosity of Sgr A* in the soft 
X-ray band is known to be variable^®. The observed fluctuations can be explained by changes in M 
by a few tens of percent over timescales of years. Variability may also arise from fluctuations in the 
parameter /3, due to random fluctuations in the magnetic field strength. The near-IR signal in the CI 
peak should show lower amplitude variations; we are not aware of any measurements. Since the S, CI 
and C2 peaks are closely coupled in the model, variations in their heights should be correlated in time, 
a testable prediction. 

The bremsstrahlung peak B occurs at the thermal frequency, i> = kTg/h, so its position is determined 
primarily by Tg. Figure 1 suggests that the model has this peak too close to the hard X-ray limits. 
With a higher Tg the peak will move to a higher i> and the discrepancy can be removed. Gamma-ray 
observations of Sgr A* at energies ~ 1 MeV will be able to constrain the properties of the bremsstrahlung 
emission and through this Tg. The height of the B peak varies as M^ since bremsstrahlung, being a 
two-particle process, is proportional to p^ . 

We briefiy mention the effect of the parameters a and /3. For a < 1, we find that models with 
a given value of M/a have virtually identical spectra. Therefore, in Fig. 1, any change in a merely 
requires rescaling M by the same factor. The effect of /3 is more interesting. As we increase /3, i.e. as 
we increase the importance of gas pressure relative to magnetic pressure, we find that Tg goes up. The 
dashed lines in Figs. 1 and 2 show an example with /3 = 0.95. This model is perhaps a better fit to 
the data than the /3 = 0.5 model. In this case, Tg max = 1-4 x 10^*^ K and the source size at 86 GHz is 
0.9 AU. 

We should emphasize that there is very little fine tuning in the models. Merely by fitting the S 
peak we are able to estimate M and M , and with these parameters we obtain a more or less satisfactory 
fit to the rest of the spectrum. The dip between S and CI is a robust prediction of the model which 
nicely explains a feature in the data very difficult to fit with other models. Further, we must emphasize 
that the model does not invoke a number of different radiating regions. All the emission comes from 
the same hot electrons (over a range of R), and the different peaks in the spectrum are just signatures 
of different radiation mechanisms. Considering how simple the model is in its basic structure, the level 
of agreement with the data is rather good, indicating perhaps that the model is correct in its essentials. 

The accreting material reaches relativistic temperatures in our model because it falls into a rela- 
tivistic potential. The only objects known in astrophysics with such deep potentials are neutron stars 
and black holes, and of the two, neutron stars are believed to be limited to masses < 3Mq. Since the 
accreting object in our model has a mass of IO^Mq, this implies that Sgr A* must be a black hole. 

While such an argument based on mass is quite plausible and is routinely employed in searches 
for black holes, we note that the crucial feature which distinguishes a black hole from other objects is 
not mass but rather the existence of a horizon, a surface which separates the interior of the black hole 
from the rest of the universe. The evidence for a black hole would be much more secure if one could 



prove that the object in question has a horizon into which matter can fall but out of which nothing, 
not even light, can escape. Remarkably, we believe that our model of Sgr A* comes close to supplying 
such a proof. 

Recall that the most important aspect of our model is that it is advection-dominated. This is the 
reason that we are able to invoke a fairly large accretion rate, M ~ few x lO~'^MQyi~^ , and yet have a 
total luminosity as low as ~ 2 x lO^^ergs""'^, much less than Mc'"^ ~ lO'^-'^ergs""'^ . Where does the rest of 
the energy go? If Sgr A* were a normal object, then all the advected energy would finally be radiated 
from its surface. This energy will presumably come out somewhere in the electromagnetic band and will 
have a net luminosity ~ few x 10'*'^ ergs" ^. There is no evidence at all that Sgr A* is emitting this much 
radiation. In contrast, in our model of Sgr A* , something like 99.9% of the energy is advected with the 
gas and disappears through the horizon and only 0.1% is radiated (see 1 — / in Fig. 2). It is this feature 
above all else that allows the model to work. Therefore, if one accepts the model we have proposed, 
then one automatically has to accept the fact that the central object in Sgr A* actually does have a 
horizon. To our knowledge this is the first "proof" of the reality of black hole horizons. Unfortunately, 
the proof is not quite complete because, in principle, the missing energy could be emitted in the form 
of kinetic energy in an outfiow. Although there is no evidence, as far as we know, for such an outfiow, 
it is probably hard to rule one out. 

The model we have presented for Sgr A* can be applied to other black hole candidates accreting 
at low rates. The majority of ultra-soft X-ray transients (also known as X-ray novae) in the Galaxy 
are believed to be stellar-mass black holes (M ~ IOMq)^^. In outburst, these sources accrete at close 
to MEdd- Under these conditions, cooling is efficient, the accreting gas probably forms a thin accretion 
disk, and the total luminosity is expected to be ~ O.lMp . However, in quiescence, the X-ray transients 
have M <^ Msdd, and in this phase the sources are possibly in an advection-dominated state. We ought 
to be able to fit the spectra of these systems with the same model we have proposed for Sgr A* but 
with different choices of M and M . AGN in their less active mode should again be similar to Sgr A* , 
though with larger masses, M ~ 10® — IO^Mq. It has been suggested that neighbouring galaxies like 
M31 and M32 harbour supermassive black holes^'*'^®. If these black holes are accreting at all, they 
are probably doing so in an advection-domination state and they should have weak radio and sub-mm 
emission corresponding to the S peak. The CI, C2 and B peaks in the spectrum will also be present, 
but since they depend more sensitively on M they will be undetectable unless M is quite high. 

Among Galactic X-ray binaries, it has turned out to be difficult to distinguish between black hole 
and neutron star systems because at high M both have similar luminosities ~ O.lMc^. However, at 
low M , if the accretion occurs in an advection-dominated mode, the difference between an object with 
a horizon and one with a surface is quite large. In the case of a low M accreting neutron star we 
expect that, even if the accretion fiow is advection-dominated, all the gravitational energy released in 
the accretion will eventually be radiated from the stellar surface, and the total luminosity will continue 
to be ~ O.lMc^. In contrast, a black hole in a similar situation will release a far smaller fraction of 
the rest energy of the accreting material. This distinction between neutron star and black hole systems 
could perhaps be used to discover new black holes in the Galaxy. 
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Figure 1. Filled circles indicate measurements of the spectrum of Sgr A* at various frequencies. The 
arrows represent upper limits. The solid and dashed lines show the spectra corresponding to two models 
described in the paper. 
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Figure 2. The variation with radius of the ion temperature, Ti, the electron temperature, Tg, and the 
radiative efficiency, 1 — /, for the two models shown in Fig. 1. 



